arXiv:1504.01736v2 [hep-ph] 10 Aug 2015 
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For potentials with n-Higgs-boson doublets stability, electroweak symmetry breaking, and the 
stationarity equations are discussed in detail. This is done within the bilinear formalism which 
simplifies the investigation, in particular since irrelevant gauge degrees of freedom are systematically 
avoided. For the case that the potential leads to the physically relevant electroweak symmetry 
breaking the mass matrices of the physical Higgs bosons are given explicitly. 


1. INTRODUCTION 


been discussed in [2, 5, 19-24]. 


Despite the fact that the Standard Model (SM) has 
only one Higgs-boson doublet, there is no theoretical 
restiction to impose a larger number of Higgs-boson dou¬ 
blets. In particular, an extended Higgs sector opens the 
possibility of CP violation in the Higgs potential. This 
was already shown by T.D. Lee for the case of the two- 
Higgs-doublet model (THDM) [1]. 

Here we want to focus on the general n-Higgs-doublet 
potential, where we assume that all doublets carry the 
same hypercharge. The aim is to find precise conditions 
for stability, electroweak symmetry-breaking, as well as 
to give equations to find systematically all stationary 
points, in particular, to detect the global minimum. It 
was shown that this is indeed possible in the case of the 
THDM [2] as well as in the 3HDM [3]. Here we want to 
generalise these findings. We will apply the bilinear for¬ 
malism, which was developed in [2, 4] and independently 
in [5]. Let us note, that the one-to-one correspondance 
of the gauge orbits of the Higgs-boson doublets and the 
bilinears in the general nHDM was already given in [2]. 

On the experimental side there is lots of effort spent to 
detect more than one physical Higgs boson for instance 
by the current LHC experiments. On the theoretical side 
also many models have been proposed which involve an 
extended Higgs sector. It is well known that supersym¬ 
metric models like the minimal and the next-to-minimal 
supersymmetric standard model require extended Higgs 
sectors. For reviews see for instance [6] and [7, 8], respec¬ 
tively. Two-Higgs-doublet models have been reviewed in 
[9]. The completely general Higgs sector was considered 
in [10] in connection with possible CP-violating effects 
in Z-boson decays. For further works on models with 
extended Higgs sectors see for instance [11-17]. Inter¬ 
esting relations between charge breaking, CP violating, 
and the normal vacuum in multi-Higgs-doublet models 
were obtained in [18]. Let us also mention that various 
aspects of the general nHDM in terms of bilinears have 
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2. BILINEARS 

Let us now consider the tree-level Higgs poten¬ 
tial of models with n Higgs-boson doublets satisfying 
SU{2)l X U{1)y electroweak gauge symmetry. The case 
of n Higgs-boson doublets is a generalisation of the cases 
with two or three doublets which were discussed in detail 
in [2, 3]. 

We will assume that we have n> 2 doublets which all 
carry the same hypercharge y = -1-1/2 and denote the 
complex doublet fields by 

; i = l,...,n. (2.1) 

The most general SU{2)l x U{1)y gauge invariant Higgs 
potential consists solely of products of the Higgs-boson 
doublets in the form 


(p,{xy(pj{x), i,j e n}. (2.2) 

We will now introduce gauge invariant bilinears, which 
turn out to be convenient to discuss the properties of the 
Higgs potential such as its stability and its stationary 
points. 

To this end we introduce the n x 2 matrix of the Higgs- 
boson fields 


fpt Pi\ (P 





(2.3) 


All possible SU{2)l x U{1)y invariant scalar products 
may be arranged into the hermitian n x n matrix 



/p’Wi 



K = = 



■ PiP2 



pIpu ■ 
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A basis for the n x n matrices is given by the ma¬ 
trices 


with all Ki > 0. We define for any hermitian matrix if 
with eigenvalues ki ,..., the symmetric sums 


Xa, a = 0,1,... — 1 (2-5) 

where 



is the conveniently scaled unit matrix and Xa, a = 
— 1, are the generalised Gell-Mann matrices. 
An explicit construction and numbering scheme of the 
generalised Gell-Mann matrices is given in appendix A. 
We will here and in the following assume that greek in¬ 
dices {a, 13, ...) run from 0 to — 1 and latin indices 
(a, b, ...) from 1 to — 1. We find 


50 := 1, 

51 := Ki+ K2 3 --I- Kn, 

32 '■= ^ 

( 2 . 12 ) 

3k •— ^ ^ ^12 ' ' ' ^ik 1 

l<ii<i2<---<ifc<n 

s„ := Ki • K 2 • ■ • ■ • Kn = det(if). 

The hermitian matrix ^ is positive semidefinite if and 
only if 

Sfc > 0 for fc = 0,..., n. (2.13) 


tr(AaA^) — 2Sq^i3, 

tr(Aa) = (5a0- 

(2.7) 

The decomposition of if (2.4) reads now 


K = iif„A„ 

(2.8) 

where the real coefficients Ka are given by 


Ka=K*= tliKXa). 

(2.9) 

Note that in particular 


ifo = tr(if Ao) = (<i>i<i>i + ■ ■. + ■ 

(2.10) 


Suppose the matrix if has rank 0, then, clearly, all Ki 
have to vanish, corresponding to 

Si = S 2 = ... s„ = 0. (2.14) 

Vice versa, starting with the conditions (2.14) for a her¬ 
mitian matrix the last condition s„ = 0 requires that 
one eigenvalue has to vanish, for instance, «;„ = 0, with¬ 
out loss of generality. The next-to-last condition in turn 
requires that another, say Kn-i = 0, and so on. There¬ 
fore we get if = 0. 

Next suppose the hermitian matrix if has rank 1, then, 
without loss of generality, we can assume 


Kl > 0, 

K2 = . . . = Kn = 0. 


(2.15) 


With the matrix if, as defined in terms of the doublets 
in (2.4), as well as the decomposition (2.8), (2.9), we may 
immediately express the scalar products in terms of the 
bilinears. 

The matrix X (2.4) is positive semidefinite which fol¬ 
lows directly from its definition X = (fxj)^. The coeffi¬ 
cients Ka of its decomposition (2.8) are completely fixed 
given the Higgs-boson fields. 

The n X 2 matrix </> has trivially rank smaller or equal 
2, from which follows that this holds also for the ma¬ 
trix if. As was shown in detail in [2] (see their the¬ 
orem 5), any hermitian n x n matrix with rank equal 
or smaller than 2 determines the Higgs-boson fields (pi, 
i = 1,... ,n uniquely, up to a gauge transformation. 

Let us now discuss the properties of the matrix with 
respect to its rank. Since the nxn matrix if is hermitian 
and positive semidefinite we can, by a unitary transfor¬ 
mation U, diagonalise this matrix. 



( Hi 

0 • 

■ 

UKU'^ = 

0 

H2 ■ 

• 0 


U 

0 • 

■ j 


If follows immediately from (2.12) 


51 > 0, 

5 2 = ... = s„ = 0. 


(2.16) 


On the other hand, having the conditions (2.16) for a 
hermitian matrix fulfilled, employing (2.12), the last 
condition = 0 requires that at least one Ki vanishes, 
for instance = 0 without loss of generality. Then the 
next-to-last condition requires that another eigenvalue 
has to vanish, for instance k„_i = 0. That is, we have 
Kn = ■ ■ ■ = K 2 = 0. Eventually, the first condition dic¬ 
tates that Kl > 0, hence, If has rank 1 and is positive 
semidefinite. 

Suppose the hermitian matrix if has rank 2, then, 
without loss of generality, we can assume 


Kl > 0,K2 > 0 , 

K3 = . . . = Kn = 0. 

If follows immediately from (2.12) 


(2.17) 


Si > 0, S2 > 0, 

S3 = . . . = S„ = 0. 


(2.18) 
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On the other hand, having the conditions (2.18) for a 
hermitian matrix if fulfilled, employing (2.12), the con¬ 
ditions S 3 = ... = s„ = 0 require that K 3 = ... = «„ = 0, 
without loss of generality. Then the first two conditions 
of (2.18) state ki -I- K 2 > 0 and ki ■ K 2 > 0, that is, we 
have Ki > 0 and K 2 > 0. Hence, if has rank 2 and is 
positive semidefinite. 

Therefore, we have shown the following theorem. 

Theorem 1. Let ^ be a hermitian matrix. 

if has rank k with fc = 0,1, 2 and is positive semidefinite 
if and only if 


Another way to parametrise all positive semidefinite 
matrices if of rank 1 and rank 2 is as follows. 

For rank 1 the matrix if has only one eigenvalue un¬ 
equal zero, say ki > 0, K 2 = ... = Kn = 0. Let u? be a 
normalised eigenvector of ^ to ki. Then we have 

K = ^KoWw\ 

w^w = 1, (2.23) 

ifo >0, Ki = J'2^o- 


So > 0,...,Sfc > 0, 

S^c-l-l = . . . = Sn = 0. 


(2.19) For the bilinears we get from (2.23) 


We may express the symmetric sums Sk defined in 
(2.12) in terms of basis-independent traces of powers 
of K. We have a recursion relation: 

So = 1 

/c 

Sk = ^ tr(if*), k = l,...,n . 

( 2 . 20 ) 

The derivation of (2.20) is given in appendix A. Explic¬ 
itly we get for /c = 1, 2,3, 

/ Ti 

Si =tr(if) = J -ifo, 

S2=\{ir\K)-iv{Ky)) 

= \{{n-l)Kl-KaKa) (2-21) 

— — (n^cKO^/30 ^a/3) ATf^^if^, 

S3 =i {i?{K) - + 2U{Ky)) . 

With the theorem 1 and (2.20) we have expressed the 
rank properties of the matrix if in terms of its eigenval¬ 
ues, respectively, traces of powers of the matrix K_. 

Based on theorem 1, (2.20), and (2.21), we can show 
that the gauge orbits of the n Higgs-boson doublet 
fields are in one to one correspondance to the vec¬ 
tors (ifo,. •., ifn 2 _i)'^ in the n^-dimensional space IR.,j 2 
satisfying 


Ka = tr(^Ao,) = \l -^ifoiyUatu. 


(2.24) 


Clearly, for any normalised vector w from (Dn we get with 
(2.23) a positive semidefinite matrix if of rank 1. 

For rank 2 the matrix if has exactly two eigenvalues 
larger than zero. Without loss of generality we can set 


= \l '■^Kosiv?{x), 

ifo > 0, 0 < X < ^. 


(2.25) 


Let Wi and W 2 be orthonormal eigenvectors of if to k\ 
and K 2 , respectively. We have then 


K = ^ ^KQ\^sm‘‘{x)wiw\ + cos^{x)w2wl ), 

wjwj = Sij, 

ifo > 0, 0 < X < 

Here we get for the bilinears 


(2.26) 


Ka =tr(ifAa) = Y'^ifo^sin^(x),ifiAaiui 

+ cos^{x),wlXaW2^. (2.27) 


Si > 0, S2 > 0, 

S 3 = ... = = 0. 


( 2 . 22 ) 


Here the Sfc, fc = 1,..., n, are constructed from the ma¬ 
trix if = KaXaf2 according to (2.20), (2.21). That is, to 
every gauge orbit of the Higgs-boson fields corresponds 
exactly one vector (ifa) satisfying (2.22) and vice versa. 
The first two relations of (2.22) are analogous to the light 
cone conditions of the THDM; see (36) of [2]. The re¬ 
maining relations in the case n > 2 are specihc for the 
nHDM. 


Clearly, the reverse also holds. For any two orthonor¬ 
mal vectors wi and W 2 the construction (2.26) gives a 
positive semidefinite matrix of rank 2. 

With (2.23) and (2.26) we have the general parametri- 
sation of all positive semidefinite matrices of rank 1 and 
rank 2, respectively. The parametrisations of the cor¬ 
responding bilinears are given in (2.24) and (2.27), re¬ 
spectively. Based on the bilinears we shall in the follow¬ 
ing discuss the potential, basis transformations, stability, 
minimization, and electroweak symmetry breaking of the 
general nHDM. 
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3. THE nHDM POTENTIAL AND BASIS 
TRANSFORMATIONS 

We now write the nHDM potential in terms of the 
bilinear coefficients, Kq, Ka, a = 1,..., — 1, 

V = ^oKo+^aKa + riQQKQ -\-2KQriaKa+ KaTjabKh, (3.1) 

where the n^{n^ + 3)/2 parameters ^ 0 i Cai Vooi Va and 
Vab = Vba are real. The potential (3.1) consists of all 
possible linear and quadratic terms of the bilinears, cor¬ 
responding to quadratic and quartic terms of the Higgs- 
boson doublets. Terms of higher order should not appear 
in the potential with view of renormalizability. Moreover, 
any constant term in the potential can be dropped and 
therefore (3.1) is the most general nHDM potential. We 
introduce the notation 

K = (Ai,..., ^ = (6, • • •, 

'n (^1: ■ • ■ 7 — 7 ^ 2) 

fp \ _ fvoo Vb \ _ fm 

IM-U- W"U e )- 

With this we can write the potential (3.1) as 

V = io.K^ + K^E^pKp. (3.3) 

Now we consider a change of basis of the Higgs-boson 
fields, <~pi{x) —?7 with 


Under the replacement (3.6), the Higgs potential (3.1) 
remains unchanged if we perform an appropriate simul¬ 
taneous transformation of the parameters 

^o=fo 7 ^' = R{U)t 

Voo=mo, ri' = R{U)ri, (3.9) 

E' = R{U)ER^{U). 

A realistic n-Higgs-doublet model contains besides the 
Higgs potential kinetic terms for the Higgs-boson dou¬ 
blets as well as Yukawa couplings which couple the Higgs- 
boson doublets to fermions. 

Under a basis transformation, that is, a transforma¬ 
tion (3.4) of the Higgs-boson doublets, or in terms of the 
bilinears, a transformation (3.6), the kinetic terms of the 
Higgs doublets are kept invariant. But in general the 
Yukawa couplings are not invariant under such a change 
of basis. 


4. STABILITY OF THE nHDM 

Now we study stability of the general nHDM poten¬ 
tial (3.1), given in terms of the bilinears Kq and K on 
the domain determined by (2.22). This is done in an 
analogous way to the cases with n = 2,3, t hat i s the 
THDM and the 3HDM; see [2, 3]. The case \fnj2Kii = 
^ ^ =0 corresponds to vanishing Higgs- 

boson fields and U = 0. For Kq > 0 we define 


Vi(a:)^ 


= U 


\'^n{x)'^ 


(3.4) 


where U £ U (n) is a n x n unitary transformation, that 
is, UW = 1„. From (3.4) we find (f>'ix) = U(j){x), and 
for the matrix A (2.4) and the bilinears 

A'(x) = UK(,x)U\ (3.5) 


fc =^7 ka = ^, fc = (fcl,...,fc„2_i)'^. (4.1) 

Now we write the rank conditions of theorem 1 in terms of 
k. With help of (2.8) we see that A = Aq- (Ao-l-A:aAa)/2. 
Therefore, the expressions Sk (2.20) are proportional to 
Kq. We define the dimensionless expressions Sk by 


Sk = 


Sk 
Rk ’ 


(4.2) 


iF'(x) = Kq{x), a'(x) = Rab{U)Kb{x). (3.6) 

Here Rab{U) is defined by 

U^XaU = Rab{U)Xb. (3.7) 

The (n^ ~ 1) x — 1) matrix R{U) has the properties 

R*{U) = R{U), R^{U) R{U) = det {R{U)) = 1, 

(3.8) 

that is, R{U) £ SOiri^ — 1). Let us note that the R{U) 
form only a subset of SO(n^ — 1). 

A pure phase transformation, U = exp(iQ:)lji, plays 
no role for the bilinears. We will, therefore, consider 
here only transformations (3.4) with U £ SU{n). In the 
transformation of the bilinears (3.6) Rab{U) is then the 
(n^ — 1) X (n^ — 1) matrix corresponding to U in the 
adjoint representation of SU(n). 


and get from ( 2 . 20 ) 


50 =17 

51 tr(Ao -H kaXa) = 

52 = 2 ^2 + kaXa) — — tr([Ao -I- fcaAa]^) ) = 


1 


(n - 1 - kaka) 7 


Sk=T^{-^y ^Sk-ztT 


Aq + kaXa 


k = 1,... ,n . 


(4.3) 


In terms of the ka we have for k the domain 'Dr. 


52 > 0 , 

5 3 = S 4 = ... = s„ = 0 . 


(4.4) 
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The domain boundary, ^'Dk^ is given by 

S 2 = ^ (n - 1 -/Cafca) = 0. (4.5) 

From (3.1) and (4.1) we obtain, for Ka> 0,V = V 2 + 
V 4 with 

F 2 = ifo J2(fe), J2(fc) := Co + (4.6) 

VA = KlJi{k), Jiik) :=r]oo + 2r]'^k + k^Ek (4.7) 

where we introduce the functions J 2 {k) and Jiik) on the 
domain (4.4). 

Stability of the potential means that it is bounded from 
below. The stability follows from the behaviour of V in 
the limit Kq —>■ 00 , hence, by the signs of J 4 (fc) and J 2 {k) 
in (4.6), (4.7). For a model to be at least marginally 
stable, the conditions 

Mk) > 0 or 

Ji{k) = 0 and J 2 (fc) > 0 

for all k G that is, all k satisfying (4.4) are necessary 
and sufficient, since this is equivalent to F > 0 for Kq — >■ 
00 in all possible allowed directions k. The more strict 
stability property F —>■ 00 for Kq —>■ 00 and any allowed 
k requires F to be stable either in the strong or the weak 
sense. For strong stability we require 

Jiik) > 0 (4.9) 


0<x<J. (4.13) 

We have to find the stationary points of 

J4{k) = kaEapkp (4.14) 

under the constraints (4.12) and (4.13). The variation 
is with respect to the real and imaginary parts of the 
components of Wi and W 2 and to %. It is easy to check 
that the gradient matrix of the four constraints (4.12) has 
rank 4. Therefore, we can use the Lagrange method and 
add these constraints with four multipliers to J 4 (4.14). 
We construct the function 

F{w\,Wi,wl,W2,X^Ui,U2,U3,U4) = 

J4{k) - Ui{w\wi - 1) - U2{wlw2 - 1) ( 4 , 15 ) 

- U3^{w\w2 + wlwi) - U4Y.{w\w2 - IcJlCl) 

where k = (ka) is to be inserted from (4.11). The equa¬ 
tions for the stationary points of J 4 are then obtained 
from 


X, Ul, U2, U 3 , U4) = 0, 


for 0 < X < 


TT 

4' 


(4.16) 


For the boundary value x = we have 


for all k G T>k] see (4.4). For stability in the weak sense 
we require for all k G T>k 


Jiik) > 0 , 

J 2 {k) >0 for all k where J 4 {k) = 0 . 


(4.10) 


In order to check that J 4 {k) is positive (semi-)definite, it 
is sufficient to consider its value for all stationary points 
on the domain Dk- This is true because the global mini¬ 
mum of the continuous function J 4 {k) is reached on the 
compact domain 'Dk, and since the global minimum is 
among the stationary points. 

In order to find the stationary points of J 4 (fc) in the 
interior of the domain Vk we note that here k of (4.1) has 
rank 2. Therefore, k and ka can be represented as shown 
in (2.26) and (2.27), respectively, but setting Kq = 1. 
This gives 


F(w\,wi,wl, W2, 7 r/ 4 , Ul, U2, U3, U4) = 0. (4.17) 

The gradients of F with respect to mi and W 2 give the 
hermitian conjugate of the gradients with respect to 
w\, respectively, in (4.16) and (4.17), thus, nothing new. 
For k, (4.1), of rank 1 we use (2.23), (2.24) to get 



where we have the constraint 


w^w — 1 = 0 . 


(4.18) 


(4.19) 


ka = 


sin^(x)m{Aami -f cos^(x)mjAam 2 (4.11) 


where 

w\wi — 1 = 0 , 

w\w2 — 1 = 0 , 

^{w\w2 + wlwi) = 0 , 
^.{w\w2 - wlwi) = 0 , 


(4.12) 


We easily check that the gradient matrix of the constraint 
has here rank 1. Therefore we add (4.19) with one La¬ 
grange multiplier to J 4 and get 

^( 11 ;^, w, u) = J 4 {k) — u{w^w — 1 ) = 

XawEafiW^ — u{w'^ w — 1). (4.20) 

The equations determining the stationary points of J 4 {k) 
on the boundary dT>k, that is, for k of rank 1 , are then 

'^■wXnF{w\w,u) = 0 . 


(4.21) 
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All stationary points obtained from (4.16), (4.17), and 
(4.21) have to fulfill the condition J^ik) > 0 for stability 
in the strong sense. If for all stationary points we have 
JAik) > 0, then for every solution k with Ji{k) = 0 we 
have to have J 2 ik) > 0 for stability in the weak sense, 
or at least J 2 {k) = 0 for marginal stability. If none of 
these conditions is fulfilled, that is, if we find at least one 
stationary direction k with Ji{k) < 0 or Ji{k) = 0 but 
J 2 {k) < 0, the potential is unstable. 

Our discussion above of the stability conditions for the 
nHDM potential generalises the results for the THDM 
and the 3HDM in [2] and [3], respectively. We have been 
careful to use in our present paper a compatible notation. 
The stability properties of the general THDM and 3HDM 
potentials were analysed in detail in chapters 4 of [2] and 
[3], respectively. Also explicit examples of THDM and 
3HDM potentials, using conventional parametrisations, 
were discussed in these references. 


5. ELECTROWEAK SYMMETRY BREAKING 
IN THE nHDM 


Now we assume that the nHDM potential is stable, 
that is, it is bounded from below. This means that the 
global minimum will be among the stationary points of V. 
We now want to distinguish the different cases of min¬ 
ima with respect to the underlying electroweak symme¬ 
try. We shall in the following present the corresponding 
stationarity equations. 

We have seen in section 2, that the space of the Higgs- 
boson doublets is determined, up to electroweak gauge 
transformations, by the space of the hermitian nx n ma¬ 
trices A with rank smaller or equal 2. Based on the fact 
that the rank of the matrix A is equal to the rank of 
the Higgs-boson held matrix cj) (2.3) we can distinguish 
the different types of minima with respect to electroweak 
symmetry breaking as follows. We start with writing at 
the global minimum, that is, the vacuum conhguration, 
the n X 2 matrix of the Higgs-boson helds as 


(<^) 


^ Vi Vi 


(5.1) 


Suppose, this matrix has rank 2, then we cannot, by 
a SU(2)l X U{ 1 )y transformation, get a form with all 
charged components u)*", i = 1,... ,n vanishing. Hence, 
the SU{2)l X U{ 1 )y group is fully broken. Next, sup¬ 
pose that at the global minimum the matrix {(p) has rank 
one. Then we can, by a SU{2)l x U{ 1 )y transformation 
get a form with all charged components u)*" vanishing. 
Further, we can identity the unbroken C/(l) gauge group 
with the electromagnetic gauge group U{l)em ■ Hence, a 
minimum with rank one corresponds to the electroweak- 
symmetry breaking SU(2)l x U{\)y U{\)em- Even¬ 

tually, suppose we get a vanishing matrix at the mini¬ 
mum, {4) = 0- This corresponds to an unbroken elec¬ 
troweak symmetry. Let us note that only a minimum 


with a partially broken electroweak symmetry is physi¬ 
cally acceptable. 

We study now the matrix corresponding to 

(</>) (5.1) 


= {4){4)^ = (5.2) 

For an acceptable vacuum (</>), K „ must have rank 1. 
From theorem 1 we see that has rank 1 and is positive 
semidefinite if and only if 


fn 

= V ^ (5.3) 

(S 2 ) = . . . = (s„) = 0. 

We can bring the vacuum value {4) of rank 1, by suit¬ 
able SU{2)l X U{1)y and U{n) transformations (3.4), to 
the form 


(<(>) 


/O 0 ' 

Vo uo/A 


Vo > 0. 


(5.4) 


In a realistic model uq must be the usual Higgs-boson 
vacuum expectation value. 


Vo « 246 GeV. (5.5) 

With (5.4) we find in this basis a simple form for 
respectively Kya'- 




/O ... 0 0\ 


0 ... 0 0 

\0 ■■■ 0 v^J 


— 2 5 


(5.6) 


(Ky^) = 




(I,0,...,0,-V^^)^. 


We note that another possible choice for the vacuum 
expectation value, achievable by suitable transformations 
from SU{2)l x U{1)y and U{n) (3.4), is 




/O vo/V^ 

• ■ 5 

Vo 0 J 


Vo > 0. 


Here we get 


(vl 0 ••• 0 \ 
1 0 0 ••• 0 


Vo 0 ••• 0/ 


(5.7) 


(5.8) 


In the cases where { 4 ) of (5.1) has rank 2 or rank 0 also 
the matrix (5-2), has rank 2 or zero, respectively. 
The corresponding conditions for K „ are given explicitly 
in theorem 1 if we replace all expressions by their vacuum 
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expectation values, that is, by Kq, by K^a. and 
Si by {si). We summarise our findings for the vacuum 
expectation values to a given potential V as follows. 

Suppose {(j)) is the vacuum expectation value of the 
Higgs-boson field matrix to a given, stable, potential V 
and = Kya^al‘2- The gauge symmetry 

SU{2)l X U{\)y is fully broken by the vacuum if and 
only if 

>0, (n - l)Kl^ - K^aK.a > 0. (5.9) 

We have the breaking SU{2)l x U{1)y —t U{l)em if and 
only if 


The gradient matrix of the constraint has rank 1 as re¬ 
quired and we get the equations 

^^ "a) 0 , ^ 

Ko>0. ^ ’ 

All stationarity matrices K = KaXal‘2 of rank 2 are 
obtained from the following system of equations. We rep¬ 
resent I£ of rank 2 and the corresponding Ka as in (2.26) 
and (2.27), respectively, and take the constraints for Wi 
into account with the help of four Lagrange multipliers; 
cf. (4.12), (4.15). We have then to determine the sta¬ 
tionary points of the function 


K,o>0, {n-l)K^o-K^aK,a=0. (5.10) 

We have no breaking of SU{2)l x t/(l)y if and only if 

= 0. (5.11) 

Clearly, we have always 

( 53 ) =••• = (Sn) =0. (5.12) 

6. STATIONARY POINTS 


F{w\,wi,wl, W2, X, Ko, Ui,U2, M3, Ui) = 

V{Ka) - Ui{w\wi - 1) - U 2 {wlw 2 - 1) (@. 5 ) 

- U3l-{w\w2 + wlwi) - U4-^{w\w2 - 
z 2i 

The stationarity equations are then 

VttR- 

w[,'W^,X,Ko,Ui,U2,U3,U4 

F{w\,Wi,wI,W2,X,Kq,Ui,U2,U3,U4) = 0 , 

0 < X < J, Ao > 0 . ( 6 . 6 ) 


Now suppose we have a stable potential. We shall 
study the stationarity equations with view on the elec- 
troweak symmetry breaking behavior. If the potential 
is stable, the global minimum is among the stationary 
points of V. We classify the stationary points by the 
rank of the stationarity matrix A. We will apply the 
conditions for A having rank 0 , 1 , 2 as given in theo¬ 
rem 1 and (2.23) to (2.27). 

Rank 0, that is, K_ = 0, respectively Ka = 0, a = 
0,..., — 1, corresponds to a stationary point of V with 

value V{Ka) = 0 . 

All stationarity matrices A = KaXaj^ of rank 1 are 
obtained from the following system of equations. We 
represent A of rank 1 according to (2.23). Then Aq, is 
given by (2.24) and V (3.3) by 


= ^o.Ka + KaEapKp 



where 


Ao>0, mi'll;-1 = 0. (6.2) 

Taking the constraint equation in (6.2) into account with 
a Lagrange multiplier u we get the following function to 
determine the stationary points of V with A of rank 1 


For X = TJ'/d we get 


,^^2 ,'*^3 ,^4 

A(m|,mi,mJ,m2,7r/4, Ao,ui,M2,U3,M4) = 0, 

Ao > 0. (6.7) 


The stationarity matrix A = KaXal‘2 with the lowest 
value of R(Ao,..., K^^-i) gives the global-minimum ma¬ 
trix A ^. of the potential. In general there may be degen¬ 
erate global minima with the same potential value. It was 
shown that systems of equations of the type (6.4), ( 6 . 6 ), 
and (6.7) can be solved via the Groebner-basis approach 
or homotopy continuation; see for instance [25, 26]. 


7. THE POTENTIAL AFTER SYMMETRY 
BREAKING 


In this section we present the calculation of the physi¬ 
cal Higgs-boson masses in the nHDM. Suppose that the 
potential is stable and leads to the desired electroweak 
symmetry breaking, that is, A^ has rank 1. From the 
previous discussion follows that the global minimum has 
then to be obtained from a solution of the set of equa¬ 
tions (6.4). 

Using (6.1) we can write (6.4) explicitly as follows 



+ ‘2-Eap 



Koiw'^Xpw) 


XaW 


E{w\w, Ko,u) = V{Ka) — u{w^w — 1 ). ( 6 . 3 ) 


uw = 0 , ( 7 . 1 ) 




w^w — 1=0, 


(7.2) Therefore, (7.10) gives 


+ 2Eal3^Ko{w^ Xi3w) 


(w^XaW) = 0, (7.3) 


Ko > 0. (7.4) 

Multiplying (7.1) with from left and nsing (7.2) and 
(7.3) we find 

u = 0. (7.5) 


(C(ti- 1)2 ~ *C(n-l)2 + l)ei + . . . + 

(Cn ^—3 ^Cn^—2)^71 — 1 

+ C-V^e„ = 0. (7.12) 

Since all (p are real we get as a resnlt of the stationarity 
eqnations for of rank 1 in onr basis 

Cp = o 

for p= (n- 1)^, (n - 1)^ + 1,..., (n^ - 2), -. (7.13) 


The vacuum value K ^. is solution of this system of 
equations. In the following we will always work in a basis 
where (</)) and have the forms (5.4) and (5.6), respec¬ 
tively. Furthermore, it is convenient to use instead of 
a = 0,1,..., — 2, — 1 the basis -k, 1,..., — 2, —; 

see appendix A. Thus, all indices p, a, ... run over this 
latter index set in the following. From (5.6) we find 


l£.v ~ — 


n/t— t 


2^2 

{K^p) = ( 0 , •••, 0, 


(7.6) 


= 




Un. 


Here and in the following e;, I £ {1,..., n}, are the usual 
Cartesian unit vectors in (Dn. We get now that for the 
solution vector w in (7.1) we have 


w = e 


n 


and that 


We define 



Cp — ^p 2Epa-Ky(j 

= ^p + 2Ep_Ky_ 


(7.7) 


(7.8) 


(7.9) 


With this we can write (7.1), using Kq > 0, in the basis 
-k, 1,..., — 2, — as 


Now we turn to the Higgs-held matrix (j). As stated 
above we work in the basis where, in the unitary gauge, 
the vacuum-expectation value {(p) has the form (5.4). For 
the original Higgs helds expressed in terms of the physical 
helds we get then 


p„(x) = ^(^^^°o(a:))’ 


* = 1,..., n — 1 


(7.14) 


with vq real and positive, neutral helds ho{x), H^{x), 
A°{x), as well as the complex charged helds H^{x) with 
i = 1,... ,n — 1. The negatively charged Higgs-boson 
helds are dehned by H~{x) = {H^{xp)\ Hence, we get 
in the nHDM the physical helds 


2n — 1 neutral helds: E[°{x),A'^{x),ho{x), 
n — 1 charged helds: (x). 


with * = 1,..., n — 1. It is clear that the n origi¬ 
nal complex doublets of any nHDM, corresponding to 4n 
real degrees of freedom, yield 2n — I real helds and n — I 
complex helds, with the 3 remaining degrees of freedom 
absorbed via the mechanism of electroweak symmetry 
breaking. Expressing the bilinears in the parametriza- 
tion (7.14) via (2.4) and (2.8) we can write the potential 
in terms of the physical helds (7.15). We start by expand¬ 
ing all quantities in powers of the physical helds. This 
gives for the held matrix 


CpApe„ = 0. (7.10) 

From the explicit construction and numbering scheme of 
the matrices Xp in appendix A we see that we have 

XpSn = 0, for p = +, 1,..., (n - 1)^ - 1, 
A(^_i)2e^ = ei, 

A(n—+ 

: (7.11) 

— ^n—li 

An^— 

A—Gtj, — 


Pix) = {(p) +4)^'^\x), 

^ Ht{x) -^{H^{x)+iA\{x)) ^ 

(x') = ' 

^ (^^n-i(a:)-ktA0_i(a:)) 

V 0 y 

(7.16) 

For iF(x) we get 

K{x) =K^+ (x) -k (x), 

^'^\x) = (^^^yx)(()^^^’^(x). 


(7.17) 
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Explicitly we get for (x) 

= ti ^((/>}'^X+(/>^^^x) + h.c.'j 

= 0 , 


Similarly we show that the only non-zero components 
of Kp^\x) are as follows: 

for p = (n — 1)^ + 21 — 2, 

Kf\x)=-voA°i{x), (7.19) 

for p = {n — lY + 21 — 1, 

kY\x) =VoV2ho{x) 

where Z = l,...,n—1. 

For the potential we write 

V = -h + V’A) + e(4) (7.20) 

and with (3.3) and (7.17) we get 

=K^P^P + K^pEp^K^a, 

=kI,^\xYp + 2K^Y^{x)Ep,K,^, 

=KY\x)ip + 2KY\x)Ep,K,, 

+ KY\x)Ep,kY){x), 
y(3) =2KY\x)Ep,kY\x), 

=Kf\x)Ep.KY\x). 

We shall now simplify the expressions for VA\ and 
using (7.6), (7.13), and (7.19). Writing as 


y(o) 


t^Eyp^p + ^Kyp\Y,p -f 2Ep(;Kyy^'^ 

2^yp{^p + Cp) 


(7.22) 


we find with (7.6) and (7.13) 


= \Kypip = i/W-C- 

- i 2 1 


(^0 - y^^e„2_i) (7.23) 


which is the potential value at the vacuum. Next we 
consider VA). With (7.9), (7.13), and (7.19) we get 


VA) =kY\x) 

=kY\x)Cp = Q 


^p + 2,EpaKy 


(7.24) 


since for each term in the above sum either kY\x) = 0 
or Cp = 0, p = -f, 1,..., — 2, —. This result must come 


out since we are expanding around the true minimum of 
the potential. 

Finally we consider Using again (7.9), (7.13), 

(7.17), and (7.21) we can write this as 

F(2) =if(2)(x)Cp +if(i)(x)F,.if(i)(x) 

= tr (^())(^)'l'(a:)CpAp(^(^)(a;)^ (7.25) 

+ kY\x)Ep,kYYx). 

Here in the hrst term the sum runs only over p = 
-f, 1,..., (n — 1)^ — 1 due to (7.13), in the second term 
only over p,a = (n — 1)^,..., — 2, — due to (7.18) and 

(7.19). The evaluation of these terms is straightforward. 
We define the fields 

n+{x) ={H+{x),...,HY_i{x))'^ , 

n°{x)={H<^ix),...,H°_Yx)f, (7.26) 

A°{x) =(H?(x),...,H°_i(a:))^. 

Furthermore, we define (n — 1) x (n — 1) matrices 




c+^-i-+ ^2 

p=i 


restricted to ’ 

the first n — 1 dimensions 
(7.27) 


£hH = (2tlQF'(„_i)2^2Z-2,(n-l)2+2Z'-2)i 
£hA = ( —2tlQF(„_i)2_|_2;_2,(„_l)2+2Z'-l)i 
£aH = (~2?;Qii'(„_i)2_|_2;_i_(„_i)2^2i'-2)i 
£aA = (2Wo^(„-l)2+2i-l.(r!,-l)2+2i'-l) 
where 

1,1’ G {l,...,n- 1}. 


(7.28) 


We also need the (n — 1) x 1 matrices 

£h- = ( 2 A/ 2 r;QF(„_i) 2 + 2 ;_ 2 ,-) 

= (~2^(„_i)2_|_2/_2), 

^A- = ( — 2 '\/ 2 Wo-^(n_i) 2 + 2 /_i,_) 
= (2C(n-l)2-|-2i-l), 

the 1 X (n — 1) matrices 

£-H = ( 2 '/ 2 UgF_^(„_l )2 + 2 p_ 2 ) 

= (~2^(„_i)2+2i'-2); 

£-a = (—2'\/2?;o^_ („_i)2_|_2//_i) 
= (2^(n-l)2+2P-l) 

and the scalar 

£— = 4,v^E— = -2y/2^_ 

/p-— (^0 1^^2_]^^ 

+ In 


(7.29) 


(7.30) 


(7.31) 
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In (7.29) to (7.31) we used (7.13) and (7.23). With all 
this we obtain for 

1/(2) =-H+t(a;) 

1 

+ ^“"^(x), ho{x)) -J^l A°{x), 

2 \ ho{x) ) 

(7.32) 

where the mass matrix squared of the charged fields, 
is given in (7.27) and that of the neutral fields, 
is given by 

+ £hh - Im (.^^\) + £ha £h^\ 
•^n ~ ( “b ^AH Re(.^^^) + £aA £a- | . 

V £-h £-a £-) 

(7.33) 

Since we have assumed that we are dealing with the 
true vacuum, y(*^) (7.23) must be below or at most equal 
to the potential value at ^ = 0. That is, we must have 

y(°) < 0 (7.34) 

which implies, from (7.23), 

Co - < 0. (7.35) 

Usually the true vacuum is required to be below the value 
U = 0 corresponding to the trivial stationary point Ka = 
0 and then the strict inequalities must hold in (7.34) and 
(7.35). For the true vacuum the squared mass matrices of 
the physical Higgs bosons must be positive semidefinite: 

> 0, > 0. (7.36) 

Looking at we see that in general it will not lead 
to mass degeneracy of all charged physical Higgs bosons. 
This is confirmed by the study of simple examples [3]. 
For the case that we have 


£h- = fA- = 0 (7.37) 

the field h^^x) is a mass eigenfield with mass squared 
value, see (7.31), 

ml^=£__ = -^V^^). (7.38) 

In this case the field ho^x) is called aligned with the vac¬ 
uum expectation value. 


the stability of any nHDM. In case of stability, the equa¬ 
tions to determine the stationary points of the potential 
have been presented. 

Of course, there are three types of vacuum solutions: 
with complete breaking, with no breaking, and with par¬ 
tial breaking of SU{2)l x U{1)y ■ For the latter case 
- the only one of physical interest - we have investi¬ 
gated the potential after symmetry breaking. The mass 
squared of the physical Higgs bosons have been given ex¬ 
plicitly. For all these investigations we have found the 
use of the gauge-invariant bilinears very convenient. For 
numerical investigations of the stability and stationarity 
equations one has to solve polynomial equations. For 
this there are approaches available, like Groebner-bases 
or homotopy continuation, which are capable to solve 
these sets of equations. We have found that the degree 
of these polynomial equations is independent of the num¬ 
ber n of Higgs bosons; see (4.15) to (4.17), (4.20), (4.21), 
(6.3) to (6.7), and (7.1) to (7.3). But the number of 
variables increases, in essence proportional to n. To con¬ 
clude: we find it remarkable that, using the method of 
bilinears, one can get a rather good overview of the prop¬ 
erties of the potentials of the nHDM, even if at first sight 
these potentials seem to be extremely involved. 
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Appendix A: Generalised Gell-Mann matrices and 
basis transformations 


Firstly, let us present a construction of the gener¬ 
alised Gell-Mann matrices Xa of dimension n, that is, 
a = l,...,n^ — 1. We start with defining the n x n ma¬ 
trix e^ej, with a 1 in the jth row and fcth column and 0 
elsewhere. Here ej, j = 1,..., n, are the Cartesian unit 
vectors in C„ 


ei = (1, 0, ..., 0)^, 

e„= (0, ..., 0, 1)^, 


(Al) 


8. CONCLUSION 

The n-Higgs-doublet model has been studied as a gen¬ 
eralization of the THDM and the 3HDM. Stability, elec- 
troweak symmetry breaking, and the stationary points of 
the Higgs potential have been discussed. We have pre¬ 
sented explicit sets of equations allowing to determine 


In terms of these matrices we construct — 1 hermitian 
traceless matrices Xa,a = 1,... ,n^ — 1 as follows. With 
k = 1, ... ,n — 1 and j = 1,..., fc we set 

Aa = + efe+ie], for a = k'^ + 2j - 2, (A2) 

Aa =-I-iefe+ie], for a = -f 2j - 1. (A3) 
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15 
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(n-1)2+7 












r4 - 1 


FIG. 1: Numbering scheme for the generalised Gell-Mann 
matrices \a (a = 0,. .., — 1). 


orthogonal n x n matrix 


/ 5'-i-o 0 

S = {Spc) = 0 1„2_2 0 

\5_o 0 5_,„2_iy 

(\f^ 0 rA \ 

Y n y n 

0 1„2_2 0 

\ S “ -4^1 

We have with a, (3 G {0,..., — 1}, 

p,a -2,-} 

S'^S' = S'S'^= 1„2, 5 = 5'^, 

Xp — Spa Xa ; 

‘X‘3 pfj • 


(A7) 


(A8) 


In addition we construct n — 1 diagonal matrices 


(A4) 

1 < 1 < n- 1. 

Eventually, we define the matrix Aq, proportional to the 
unit matrix, 


A, 




1(1 + 1 ) 



Let us note that the matrices Aq, (a = 0,..., — 1) de¬ 

fined in this way in particular fulfill the conditions (2.7). 
An easy way to remember this numbering scheme is as 
follows. We draw an n x n square lattice and insert the 
numbers a = 0,1,..., — 1 as shown in Fig. 1. If a is 

the upper (lower) number in an off-diagonal square then 
Aq gets a 1 (—i) in this place, 1 (+i) in the transposed 
place, and zero elsewhere. If a is in a diagonal square Aq 
is given by (A4) for a > 0 and by (A5) for a = 0. 

For n = 3 the matrices Aa (a = 1,..., 8) as constructed 
above are the standard Gell-Mann matrices; see for in¬ 
stance [27]. 

Returning to the case of general n we find it convenient 
to define also 


With the help of S we transform also ATq, ^q, Eap (see 
(2.9), (3.2), (3.3)) to the basis p, cr G {+, 1,..., — 2, —} 

Kp = SpaKa = tr(KXp), 

4,p = Spa4,ai (-^ 9 ) 

Epa — SpaEapS^^. 

This gives, for instance, with (2.4) 


K+ — \ - Ko + \ —Kn'^-i 

\ n \ n 

^ + ■ ■ • + , (AlO) 

K- = \ -Kq - J - Kn2_i = 

\ n \ n 


Appendix B: Symmetric sums 


Here we want prove the recursive relation (2.20) for 
the symmetric sums as originally defined in (2.12). Con¬ 
sider 1 < fc < n. First we note that Sk(Ki ,..., k„) is a 
homogenous function of degree fc in ki ,..., Therefore 
we have 


Ki—Sk{m, ...,Kn) = kskif^l, - . . , /^n)- (Bl) 




On the other hand we have 



(A6) 


The change form the basis Aq, Ai,..., A^ 2 _ 2 , A^ 2 _i to 
A+, Ai,..., A^ 2 _ 2 , A_ is made with help of the following 


Sk 7 ■ • ■ ; ^n) 

OKI 


E 


■■■■ ^+-1 


= Sfe_l(Kl, . . . ,K„) - ( ^ K^^■ ...■ 

1 < i 1 < ... < i _ 2 < n 

= Sfe-l(Kl, ■•■,««)- Sfc_2(Kl, . . . , Kn)«^; 

+ ... + (-l)'=-iso«;f-\ 


(B2) 

































which proves (2.20). 
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Multiplying in (B2) with k;, summing over Z, and us¬ 
ing (Bl), we get 

kSk{Ki, ...,Kn)= Sfc_i(Ki, . . . + . . . + K^) 

~ Sfe_2(Ki, . . . , K^) 

(—1)^” ^so(ki -I- .. ■ -I- (B3) 
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